Resonant Adiabatic Passage with Three Qubits 
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We investigate the non-adiabatic implementation of an adiabatic quantum teleportation protocol, 
finding that perfect fidelity can be achieved through resonance. We clarify the physical mecha- 
nisms of teleportation, for three qubits, by mapping their dynamics onto two parallel and mutually- 
coherent adiabatic passage channels. By transforming into the adiabatic frame, we explain the 
resonance by analogy with the magnetic resonance of a spin-1/2 particle. Our results establish a 
fast and robust method for transferring quantum states, and suggest an alternative route toward 
high precision quantum gates. 
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Fault-tolerant quantum computation requires high- 
precision quantum gates with noise thresholds be- 
tween 10 -4 and 10 -2 , depending on the fault-tolerance 
scheme 0, Q . This stringent requirement poses signifi- 
cant technical challenges, even for the more mature qubit 
architectures, such as those based on trapped ions Q. 
Identifying gate protocols that are both fast and robust 
is therefore an important research objective for quantum 
information processing. 

One potential approach to robust quantum gates is 
based on the adiabatic principle - a fundamental tenet 
of quantum mechanics [J]. According to the adiabatic 
theorem, a quantum system in an cigenstatc remains 
there, provided that the Hamiltonian varies slowly in 
time. Applications of the adiabatic theorem include 
the Born-Oppcnhcimcr approximation and the Landau- 
Zcner-Stiickclbcrg transition at an avoided crossing, the 
latter having been demonstrated in both superconduct- 
ing and spin qubits 0, Other experimental imple- 
mentations include adiabatic population transfers be- 
tween two or three-level systems, known as adiabatic 
passages (AP) which have been demonstrated in 

atomic, molecular, and optical devices. There are also 
theoretical proposals for realizing AP with superconduct- 
ing qubits (llj and quantum dot arrays (l2j . Adiabatic 
quantum information processing [ilj ] entails the adiabatic 
transformation of the ground state of an initial Hamil- 
tonian into that of a target Hamiltonian. Compared 
to the quantum circuit model, adiabatic gates are resis- 
tant to decoherence when a finite excitation gap persists 
throughout the evolution, and they are robust to gating 
errors, by virtue of adiabaticity. This can be a drawback 
however, since the maximum speed of an adiabatic gate 
is also proportional to the spectral gap. 

In this Letter we investigate a non-adiabatic form of 
adiabatic quantum teleportation (AQT). Conventional 
AQT was proposed in the context of fault tolerant quan- 
tum computation (3). Here, we focus on systems with 
three qubits, where we can solve the evolution analyti- 
cally. We show that resonances occur, enabling telepor- 



tation that is fast, fault-tolerant, and potentially perfect. 
This surprising effect can be explained in the language of 
spin resonance, by transforming into the adiabatic frame. 
Our results point toward a new paradigm for quantum 
algorithms, based on fast adiabatic gates. Our work also 
provides an interesting mapping between three coupled 
qubits and a three-level atom, which could lead to fur- 
ther spin analogies from atomic A-systcm physics. The 
experimental requirements for implementing AQT have 
already been demonstrated in the laboratory for triple 
quantum dots Q and superconducting circuits fl5j . Our 
results could therefore be tested immediately. 
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Figure 1: (color online). The couplings and the informa- 
tion distribution (by color and arrow) among three qubits 
are shown at the initial, intermediate, and final stages of adi- 
abatic quantum teleportation. 

The adiabatic quantum teleportation protocol is illus- 
trated in Fig. [U Initially, qubit 1 is isolated and pre- 
pared in an arbitrary superposed state, while qubits 2 
and 3 are coupled, as described below, and prepared in 
the maximally entangled singlet state. The antiferromag- 
netic coupling between qubits 1 and 2 (2 and 3) is then 
turned on (off) slowly. When the evolution is complete, 
the quantum state of qubit 1 will be teleported to qubit 
3. As proposed in Ref. the scheme succeeds when 
the run time T satisfies the adiabatic theorem. 

To explore non-adiabatic effects, we solve the exact 
dynamics of the three-qubit system. It is governed by 
a time-dependent Hamiltonian, which smoothly changes 
from the initial Hamiltonian H; at t = to the final 
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Hamiltonian Hf at t = T: 

H(t) = f(t)H i + g(t)H f . (1) 

The initial and final Hamiltonians are given by 

Hi = J { a 2x &3x + <?2y &3y + 7 a 2z o Zz ) , (2a) 
H f = J ( a lx a 2x + <Ti y a 2y + 7 a lz a 2z ) , (2b) 

where <Zj„ are the Pauli operators with i = 1,2,3 and /x = 
x, y, z, and J is the strength of the qubit-qubit coupling. 
The anisotropy parameter 7 = corresponds to an XX 
coupling, available in superconducting qubits, while 7 = 
1 corresponds to the isotropic Heisenberg coupling of spin 
qubits. The interpolation or switching functions f(t) and 
git) satisfy /(0) = g(T) = 1 and f(T) = g(0) = 0. Here 
we consider two ways to connect Hi to Hf. (i) a linear 
interpolation with f{t) = 1 — t/T and git) = t/T, and 
(ii) a harmonic interpolation with /(t) = cos(7rt/2T) and 
git) = sm(7rf/2T). 

AQT begins with an initial three-qubit state given by 
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where a|0) + is the arbitrary state to be teleported. 
The success of AQT is measured by the fidelity -F(T) = 
\(ipT\ipiT))\ 2 where \tp(T)) is the final state at time T and 
W) = 75 (|01) - |10)) 12 g> (o|0) + 6|1)) 3 is the target 
state. The dynamics of AQT is governed by the time- 
dependent Hamiltonian (0, with the initial state (0. 

Hamiltonian ([1]) satisfies the commutation relation 
[H(t), S z ] = 0, so that the z-component of the total spin 
angular momentum S z = \ ip\ z + a 2z + a 3z ) is a good 
quantum number, which is conserved during evolution. 
The three-qubit Hamiltonian (JXJ) is thus block-diagonal: 



Hit) = H^®H^®Hi © H x 

up down 



(4) 



The two H\ operators act on |000) and [111), 
while the two H3 operators act on the distinct sub- 
spaces H 1/2 = Span(|100),|010),|001)) and H-x/2 = 
Span(|011), 1 101) , 1 1 10) ) , and have the same form 



H 3 (t) = J 



if - .9)7 2 5 
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(5) 



H3 is also the AP Hamiltonian for a 3- 
10) . with the switching functions fit) and 



Interesting! 
level atom 

git) being the Stokes and pump pulses in the context of 
AP. For the initial states we consider, the Hi operators 
are never involved in the system evolution. The AQT 
protocol therefore consists of two parallel, identical and 
mutually-coherent APs governed by -H3, corresponding 
to the S z = ± h components of the three-qubit system. 



To understand the dynamics of AQT, we solve the 
time-dependent Schrddinger equation with Hamiltonian 
([5]) in two ways. First, we consider the adiabatic limit, 
for which there is a mapping between AQT and two 
mutually-coherent APs. Second, we obtain numerical so- 
lutions (and in one case, an analytical solution) for finite 
T. We also consider the separate cases of XX and Heisen- 
berg couplings. 

The adiabatic theorem states that, starting from 
an eigenstate |i?„(0)), the adiabatically evolved state 
\ip{t)) ~ e-Wo E n (t')dt'+ 7B \E n (t)) is simply an instan- 
taneous eigenstate, up to a phase factor. Here, jb is 
the Berry phase, and E n and \E n ) are the instantaneous 
eigenvalues and eigenstates of H^it), defined by 



H 3 (t)\E n (t))=E n (t)\E n (t)). 



(G) 



Solving Eq. © for an XX coupling (7 = 0), gives 
the instantaneous energy levels Eo(t) = and E±(t) = 
±2Jy / f 2 + g 2 , with the corresponding eigenstates 



\Eo(t)) = 
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Here, the mixing angle 6 = tan -1 [g(i)//(£)] runs from 
to 7r/2 as time t goes from to T. For the Heisen- 
berg coupling (7 = 1), the instantaneous energy lev- 
els are E jt)/J = (/ + g) and E±(t)/J = [-/ - g ± 
2\J f 2 — fg + g 2 ], with the corresponding eigenstates 



\E (t)) = ^=[l 1 1] T 
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where Af = 2(2 cos 9 — sin 6)Jq + 4q, q = 1 — cos 9 sin 9. 

Equation (|5|) governs the evolution of both AQT and 
conventional AP, as depicted in Fig. [5] For AP, the pop- 
ulation of a A- type system is transferred from state |1) to 
1 3), while state |2) remains unpopulated fsTlloj. Paradox- 
ically, the AP pulse sequence appears to occur in reverse 
order (S followed by P), as shown in panel (a). The in- 
stantaneous eigenstate used in this evolution is \Eo{t)) 
from Eq. (0. For AQT, on the other hand, the instanta- 
neous eigenstate used is \E-) in Eqs. (0 or (JSJ, leading 
to slight differences between panels (a) and (b) . In panels 
(b) and (c), we see that the "up" state (|0)) is transferred 
from the left-most qubit to the right-most qubit, follow- 
ing a similarly counter-intuitive pulse sequence. Since 
the i?3 operators are identical for the subspaces Hi/ 2 
and H_i/2) their separate evolutions are also identical. 
Thus, as illustrated in Fig.HJ an arbitrary state a|0)+6|l) 
of qubit 1 in Eq. (0 is transmitted to qubit 3 via two 
mutually-coherent evolutions. 
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Figure 2: (color online). Adiabatic passage protocols for (a) 
A- type levels, (b) the "up" component (red |0)) of a 3-qubit 
system, and (c) the corresponding spin configurations. From 
right to left: the initial, intermediate, and final stages of evo- 
lution. In (a) and (b), filled circles represent populations of 
levels. S and P stand for the Stokes and pump pulses. 



The adiabatic solution described above is only valid 
when JT/h > 1. In this limit, Eqs. © and © give a 
perfect (adiabatic) fidelity, F^iT) = 1. When T is finite 
however, the adiabatic theorem predicts that 1 — F oc 
( JT/h)- 2 . To obtain an infidelity 1 - F < 1(T 6 , the adi- 
abatic gate time should be JT/h ~ 0(1O 3 ), much longer 
than a conventional gate, for which JT/h ~ 0(1). Such 
slow adiabatic evolution could obviously cause problems, 
despite its intrinsic fault tolerance. However, when we 
perform a numerical integration of the time-dependent 
Schrodinger equation governed by Hamiltonian ([5]), we 
find that the infidelity 1 — F as a function of evolution 
time T is far from a smooth quadratic function. Instead, 
while it approaches the predicted upper envelope, there 
are also striking resonance features where the infidelity 
dips to zero, as shown in Fig. [3] 

The origin of the unexpected resonances in AQT fi- 
delity becomes clear when we consider the XX coupling 
with harmonic interpolation functions. For this special 
case, we can obtain an analytical solution by transform- 
ing into the adiabatic frame [l(| , as illustrated in Fig. |4] 
We define 

D{t) = A-\t)H 3 {t)A{t) , m)) = A{t)W)) , (9) 

where the column vectors of A(t) are the instantaneous 
eigenstates given by Eq. ©. The Schrodinger equation 
in the adiabatic frame takes the form 



ih§- t m)) 



D(t) -ihA' 1 ^) 



dA(t) 



dt 
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(10b) 
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Figure 3: (color online). Infidelity (1 — F) is plotted as a 
function of JT/nh for (a) XX and (b) Heisenberg couplings, 
using linear and harmonic interpolation functions (solid red 
and dotted blue lines, respectively). The black dashed line 
represents (1 — F) oc 1/(JT) 2 for comparison. 



For harmonic interpolation and XX couplings, the trans- 
formed Hamiltonian H tr becomes time-independent: 



i?tr = 2J 
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= hujQ Z + huj\ Y , 



(11a) 
(lib) 



where huo = 2 J is the absolute value of the ground state 
energy and u>i = ?t/2T is the frequency of the switch- 
ing functions, f(t) = cos(wit) and g(t) = sin(wii). The 
matrix Y', which resembles the angular momentum op- 
erator I y of a spin-1 system, is responsible for the non- 
adiabatic behavior. It has the same eigenvalues as Z, 
i.e., and ±1. The Hamiltonian in the adiabatic frame, 
iJtr = h£l{Z cos a + Y' sin a), has the eigenvalues 



e = , e± = ±hfl , with ft = \/u)q + oj 2 , (12a) 
and the corresponding eigenstates 



1 

71 



— sin a 




1 =F cos a 


iV2cosa 


, Ie ± > - \ 




sin a 




1 ± cos a 



, (12b) 



where tana = cji/wo- 

As illustrated in Fig. [U the time evolution in the adi- 
abatic frame is analogous to the rotation of a spin-1 sys- 
tem around an effective, constant magnetic field given by 
$7 = loqZ + uj{Y' , where Y' is the rotation axis associ- 
ated with matrix Y'. The state vector is initially oriented 
along Z, which corresponds to \E_(t)) in the original 
frame of Eqs. (0. In the adiabatic limit Wo > wi, the 
precession axis is Cl = Z, so the state vector does not 
precess. Thus, in the original frame, the state vector is 
given by |£L(i)) for all t. 
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Non-adiabatic evolution occurs when oj\ > 0. The 
state vector is initially aligned with Z in the adiabatic 
frame; however it precesses when Cl ^ Z . As the state 
vector deviates from Z in the adiabatic frame, it also 
deviates from the adiabatic ground state \E_{t)) in the 
original frame. After a full precession period given by 
OT = 27rn, the state vector returns to the Z direction, 
or the ideal target state \E_(T)}. The physical picture is 
analogous to the magnetic resonance of a spin- 1/2 parti- 
cle in a static magnetic field, with a small perpendicular 
ac field. In this case, the state vector precesses about a 
static magnetic field in the rotating frame 171. Il8|. 




|V(0)) = |£_(0)) 



m) + \E-(t)) 



Figure 4: (color online). Schematic representation of the time 
evolution in the adiabatic frame, given by Eq. (|13[) . The in- 
stantaneous eigenvector is represented by a thick blue arrow 
and the exact evolved state by a thick red arrow. 

As depicted in Fig. [4] the evolved state in the original 
frame is given by 



m)} = A(t)U tr (t)A- 1 (0)\i>(0)), 



(13) 



where the time evolution operator in the adiabatic frame 
is given by U tr = e ~ iHt ' t/h . The fidelity at time t = T 
can be obtained exactly: 



F(T) 



i r 



cos 2 (rJT)(l + cos 2 a) 2 + 4sin 2 (OT) cos 2 



2 cos(fJT) sin 2 a(l + cos 2 a) + sin 4 a 



(14) 



The results are indistinguishable from the numerical so- 
lution shown in Fig. |3|a). Perfect fidelity occurs at the 
resonance condition OT = 27m, which is given by 



JT/h-K 



1 

16 



n , n e N . 



(15) 



We have now identified two paths to perfect tclcporta- 
tion. The first corresponds to the asymptotic (adiabatic) 
limit on the far right-hand side of Fig. [3] The second 
occurs at any one of the resonant conditions. It is in- 
teresting that resonances only occur in certain interpola- 
tion schemes. For example, the quadratic interpolation 



and g(s) 



f(s) = 1 - s 2 and g(s) 



- s(2 — s) has resonances, while 
s 2 does not. 



Our results can be tested experimentally using cur- 
rent technology. Controllable three-qubit systems have 
been demonstrated in quantum dots @ and supercon- 
ductors [l5|. Single-shot measurements and the prepa- 
ration of singlet states are almost routine [l9[. AQT 
could therefore be implemented as follows. Qubit 1 is 
initially prepared in the "up" state, while qubits 2 and 3 
are prepared in a singlet state. After switching / and g 
according to the AQT protocol, qubit 3 is measured. Re- 
peating this experiment many times provides a fidelity 
estimate for AQT, over the evolution period T. The 
resonant peaks of the fidelity can be examined in the 
time domain by varying T. Since AQT corresponds to 
two parallel APs for the two spin components of a qubit, 
we could also explore interesting phenomena like coher- 
ent population trapping and electromagnetically induced 
transparency pot ], which have also been studied in the 
context of AP, for three-level atoms. 

So far we have only explored AQT with three qubits, 
where qubits 2 and 3 are initially in a singlet state - the 
same initial state used for conventional quantum telepor- 
tation. An interesting next step would be to study AQT 
over longer distances. Our preliminary numerical studies 
suggest that AQT could be implemented in a more gen- 
eral spin chain geometry. We leave this for future work. 

In conclusion, we have shown that adiabatic quantum 
teleportation consists of two adiabatic passages corre- 
sponding to the quantum information transfer of "up" 
and "down" components of a qubit. When this protocol 
is performed non-adiabatically, resonances occur in the 
fidelity, in analogy with magnetic spin resonance. The 
observation of resonances points toward a new paradigm 
for fast and robust adiabatic gates. Our results can 
be tested experimentally using superconducting or spin 
qubits, with currently available technologies. 
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